Periodic wave packet reconstruction in truncated tight-binding lattices 
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A class of truncated tight-binding Hermitian and non-Hermitian lattices with commensurate energy 
levels, showing periodic reconstruction of the wave packet, is presented. Examples include exact 
Bloch oscillations on a finite lattice, periodic wave packet dynamics in non-Hermtian lattices with 
a complex linear site-energy gradient, and self-imaging in lattices with commensurate energy levels 
quadratically-varying with the quantum number. 
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Introduction. Tight-binding lattice models provide 
a simple tool to investigate the transport properties of 
different physical systems, such as electronic transport 
in semiconductor superlattices or coupled quantum 
dots, photonic transport in waveguide arrays or coupled 
optical cavities, and transport of cold atoms in optical 
lattices. In certain conditions, the dynamics of a 
wave packet in the lattice turns out to be periodic. 
Such a special regime is attained, for example, in an 
infinitely-extended periodic lattice with a superimposed 
dc field, leading to the formation of an equally-spaced 
Wannicr-Stark ladder energy spectrum and to the 
appearance of Bloch oscillations (BOs) 1]. Another 
important case is that of an infinitely-extended periodic 
lattice with an applied ac field, for which quasi-energy 
band collapse leading to dynamic localization can be 
found at special ratios between amplitude and frequency 
of the ac modulation @, Such periodic dynamical 
behaviors realize wave self-imaging, i.e. a periodic 
reconstruction of the initial wave packet distribution. 
Self-imaging phenomena in tight-binding lattices have 
been observed in different physical systems (see, for 
instance, [4-8]). Even in an infinite periodic lattice 
without any external applied field periodic self-imaging 
can be observed for special periodic initial wave distri- 
butions owing to the discrete analogue of the Talbot 
effect @. Unfortunately, lattice truncation, defects and 
edge effects generally destroy the periodic wave packet 
dynamics because of incommensurate frequencies of the 
eigenstates. As a consequence, the system never revives 
fully to its initial state, and phenomena like wave packet 
collapse and revivals are generally found Recently, 
a few examples of lattice engineering that enables to 
restore exact per iodic wave packet dynamics have been 
proposed [Ill4l5| . including harmonic oscillations in 
a finite lattice with engineered hopping rates [ill EH 
and Bloch oscillations on a semi-infinite lattice with 
linearly-increasing hopping rates (14]. In this Report we 
introduce a more general class of truncated tight-binding 
lattices -not necessarily Hermitian- that sustain a 
periodic dynamics of the wave packet and that include 



the Hermitian lattices of Refc.flll Il2l . Ii"4j as special 
cases. In particular, we prove the existence of Bloch 
oscillations on engineered finite lattices, either Hermitian 
or non-Hermitian, and show self-imaging phenomena in 
a novel class of finite lattices with commensurate energy 
levels that depend quadratically on the quantum number. 

Lattice synthesis. Let us consider a tight-binding 
one-dimensional lattice with site energies V n and hop- 
ping amplitude K n between adjacent sites \n) and |n + l). 
In the nearest-neighboring approximation, indicating by 
^>(n) the occupation amplitude at lattice site \n), the 
energies E of the lattice Hamiltonian H are found as 
eigenvalues of the second-order difference equation 

%ip{n) = K n -itjj(n— 1) + K, n t/j(n+ 1) + V n ip(n) = Eip{n). 

Lattice truncation is realized by assuming K-\ — for a 
semi- infinite lattice comprising the sites |0), |1), |2) or 
K-i = km = for a finite lattice comprising the (N + 1) 
sites |0), |1), |2), ... \N). Note that the Hamiltonian H 
is Hermitian provided that V n and real-valued, 
the emergence of complex values for either hopping rates 
or site energies being the signature of a non-Hermitian 
lattice. Our aim is to judiciously engineer the hopping 
amplitudes n n and site energies V n in such a way that 
the lattice energies E form a set of commensurate num- 
bers, thus ensuring a periodic dynamics of an arbitrary 
wave packet ip(n,t) that evolves according to the time- 
dependent Schrodinger equation i(dip/dt) = T-Lip. Exam- 
ples of such a lattice engineering, corresponding to a set 
of equally-spaced discrete energy levels, have been previ- 
ously presented in Refs.[ll], 0, E3|- Here we show that 
such previous examples belong to a more general class of 
lattices with commensurate energies, which can be syn- 
thesized using the factorization method (lq|. or the so- 
called supersymmetric quantum mechanics [T7[ , adapted 
to 'discrete' quantum mechanics and difference equations 
(see, for instance, [HJ]). To this aim, let us assume that 
the hopping amplitudes k„ and site energies V n can be 
derived from two functions Fi (n) and F 2 (n) according to 
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«n - VFi(n)F 2 (n + l) , V n = - [F x (n) + F 2 (n)}. (2) 
In this case, the Hamiltonian T-L can be factorized as W = 



2 



AB, where [H] 

A = V^expQ^-v^Rexp^A) (3) 

To realize a finite lattice, with nonvanishing amplitudes 
ip(n) at the (7V+1) lattice sites |0), |1), \N), we require 
Fi(N) — and F 2 (0) = 0. For a semi-infinite lattice, 
we require F 2 (0) = solely. Let us then indicate by 
ipo(n) the solution to the difference equation Bipo(n) = 
0, i.e. y/ F 2 (n + l)i/> (n + 1) = y /F\{nj'ipo(n) ■ It then 
follows that Hipa — 0, i.e. E = is an eigenvalue of % 
[l9| . To find the other eigenvalues of T-L , let us search for 
a solution to Eq.(l) in the form ip(ri) — ipo(n)Q(n); it 
then readily follows that Q(n) satisfies the second-order 
difference equation 

F 1 (n)Q(n+l)+F 2 (n)Q(n-l) = [F 1 {n)+F 2 (n)+E}Q(n). 

(5) 

Assuming for Fi(n) and F2(n) a linear or quadratic 
functions of n, Eq.(5) can be solved in an exact way 
by assuming for Q(n) a polynomial of degree M, i.e. 

Q{ n ) = 12fLo a i nl - Substitution of such an expression 
of Q(n) into Eq.(5) yields the following relation 




M 

= E^2a p n p (6) 

p=0 

which determines the eigenvalues E and the coefficients 
a p of the polynomial Q{n) after comparison of the terms 
of the same power n p on the left- and right-hand-sides 
in Eq.(6). In particular, the eigenvalue E is found by 
comparison of the highest-order power terms n M . It is 
worth discussing in a separate way the cases of linear 
and quadratic forms for Fi(n) and F2(n), which yield 
a different dependence of the energy eigenvalues on the 
quantum number M. 

Lattices with equally- spaced energy levels: Block oscil- 
lations on a finite lattice. Let us consider a finite lat- 
tice comprising the sites n = 0,1, 2, ...,7V and assume 
Fi(n) — a(n — N) and F2(n) = 771, with a and 7 arbi- 
trary parameters with 0^7. The hopping amplitudes 
and site energies of this lattice are thus given by [see 
Eq.(2)] 

K n = \foq(n - N)(n + 1) , V n =aN - (a + -f)n. (7) 

Note that the lattice Hamiltonian is Hermitian for 07 < 
0, whereas it is not Hermitian for 0:7 > or for complex- 
values for a and/or 7. Substitution of such expressions 
of Fi(n) and F2(n) into Eq.(6) and by comparing the 
coefficient of the highest-degree term (n M ), yields the 
following simple expression for the eigenvalues 

E = M{a- 1 ), (8) 



with M = 0, 1, 2, N. Owing to the equal spacing of 
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FIG. 1: (color online) Exact Bloch oscillations in an Her- 
mitian truncated lattice, (a) coupling amplitudes Kn, and 
(b) lattice site energies V n as given by Eq.(7) for iV = 15, 
a — 1 and 7 = —0.2. (c) Numerically-computed temporal 
evolution of the site occupation probabilities \tp(n,t)\ 2 in the 
(n, i) plane for an input Gaussian wave occupation ampli- 
tudes ^(n,0) = exp{-[(n- iV/4)/4] 2 }. The inset in (c) shows 
the detailed temporal evolution of the occupation probability 
I t/j (0, t)\ 2 of the left-edge (n = 0) lattice site. 

energy eigenvalues, the lattice model (7) yields a periodic 
temporal dynamics with period T = 27r/|a — 7|. Note 
that, for a — —7 real- valued, this lattice model corre- 
sponds to the one previously introduced in Refs.fTTI. [T^|. 
i.e. to uniform site energies V n and to an optimal 
coupling of adjacent lattice sites. Interestingly, for real 
values of a and 7 but with a 7^ —7 and < 0, a linear 
gradient of site energies V n occurs, and the lattice model 
(7) thus sustains exact BOs on a finite lattice. Figure 
1 shows, as an example, a typical behavior of hopping 
amplitudes and site energies for such an Hermitian 
lattice, together with the characteristic BOs behavior for 
an initial Gaussian wave packet as obtained by numerical 
simulations of the Schrodinger equation i(dip/dt) = Htp. 
To show that exact wave packet reconstruction is insen- 
sitive to edge effects, in the example of Fig.l wave packet 
excitation was intentionally chosen close to one edge of 
the lattice, however periodic wave packet reconstruction 
is observed for any initial wave packet, regardless of its 
shape or position in the lattice. It should be also noticed 
that periodic wave packet reconstruction does not mean 
nor imply shape-invariance of the wave packet evolution. 
Remarkably, a periodic quantum evolution, associated 
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FIG. 2: (color online) Bloch oscillations in a non-Hermitian 
truncated lattice, (a) coupling amplitudes K n , and (b) lattice 
site energies V n (real and imaginary parts) as given by Eq.(7) 
for TV = 15, a = 1 + 0.2i and 7 = -1 + 0.2i. (c) Numerically- 
computed temporal evolution of the site occupation proba- 
bilities \ip(n,t)\ 2 in the (n,t) plane for an input Gaussian 
wave occupation amplitudes ^(n, 0) = exp{ — [(n — 7V/4)/4] 2 }. 
(d) Detailed temporal evolution of the occupation probability 
|^(0, t)\ 2 of the left-edge (n = 0) lattice site, (e) Behavior of 
total occupation probability P(t) = ^2 n \ip(n,t)\ 2 . 



to the existence of equally-spaced real-valued energies, 
is found even for non-Hermitian lattices by taking 
a = a + ip and 7 = — a + ip, where p and a are arbitrary 
real-valued and nonvanishing constants. In this case, 
according to Eq.(7) the hopping amplitudes K n turn 
out to be real- valued [n n = a/ (p 2 + a 2 )(N — n)(n+ 1)], 
whereas the site energies V n are complex-valued and 
linearly varying with the index n [V n = (<r + ip)N — 2ipn], 
i.e. the linear gradient of site energies is now imaginary 
rather than real-valued as in ordinary BO lattice models. 
Physically, such a new class of non-Hermitian lattices 
could be realized by a sequence of (N + 1) active optical 
waveguides with a judicious engineering of waveguide 
spacing (to tune the hopping rates K n ) and with a 
controlled gain/loss coefficient V n (see, for instance, 
[H). In spite of non-Hermiticity, the energy levels are 
real-valued and equally spaced like in the Hermitian 
lattice model. However, the kind of periodic wave packet 
evolution in the such a non-Hermitian lattice is very 
distinct from ordinary BOs found in Hermitian lattices 
(see e.g. @). In particular, oscillations of the total 
occupation probability P(t) — ^2 n \^{n,t)\ 2 are found 
as a result of the non-conservation of the norm, like in 
other non-Hermitian crystals (see for instance |21(). As 
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FIG. 3: (color online) Self-imaging in an Hermitian lattice 
with commensurate energy levels that depend quadratically 
on the quantum number [Eqs.(9) and (10)]. Parameter values 
are TV = —a — 10 and 7 = 0. (a) Behavior of k„. (b) Behavior 
of V n - (c) Numerically-computed temporal evolution of the 
site occupation probability \tp(0, t)\ 2 of the left-edge (n = 0) 
lattice site for initial single-site excitation tjj(n, 0) = <5 n , o- 



an example, Fig. 2 shows the behavior of n n and V n for 
a non-Hermitian lattice defined by Eq.(7), together with 
a typical temporal evolution of an initial Gaussian wave 
packet, showing a characteristic oscillation of the total 
occupation probability. 

As a final comment, it is worth mentioning that BOs on 
a semi-infinite Hermitian lattice with linearly-increasing 
hoping amplitudes, recently predicted in Ref . [14j . 
can be found by assuming Fi(n) — a(n + 1) and 
■^2(71) = 771, with aj > (to ensure lattice Hermiticity) 
and I a I < I7I (to ensure boundness of ?po(n) as n — > +00 
19]). In this case one obtains n n = J(n + 1) and 
V n = —fn — a (n = 0,1,2,3,...), where J = ^/cry 
and / = (a + 7) are the hopping rate and site energy 
gradients, respectively. Note that the condition |/| > 2 J 
is always satisfied, i.e. the existence of a Wannier-Stark 
energy level spectrum requires a minimum value of the 
dc force |/|. As discussed in 14], such a minimum 
value \ f\ m = 2J of the dc force corresponds to the exis- 
tence of a metal-insulator transition in this lattice model. 

Self-imaging in lattices with commensurate energy 
levels guadratically-varying with the quantum number. 
As a second class of lattices showing self-imaging phe- 
nomena, let us consider a finite lattice comprising the 
sites n = 0, 1, 2, N and assume Fi(n) — (n— N)(n+a) 
and F2(n) — n(n + 7), where a and 7 are real- valued 
parameters. The hopping amplitudes and site energies 
of this lattice are thus given by [see Eq.(2)] 

K n = V(n + a )( n ~ N )( n + + 7 + 1) (9) 
V n = -[2n 2 + (a + 7 - N)n - aN]. (10) 
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Substitution of such expressions of F\ (n) and F2 (n) into 
Eq.(6) and by comparing the coefficient of the highest 
degree term (~ n M ) in the power expansion, one obtains 
the eigenvalues 

E = M 2 -M(N + 1 + 7-a), M — 0, f , 2, N, (11) 

i.e. the energy eigenvalues are described by a quadratic 
function of the quantum number M. Provided 
that (N + f + 7 — a) is a rational number, i.e. 
(N + f + 7 — a) = T\Iti with r% and r-i integer numbers, 
the energy levels are commensurate numbers and the 
wave packet dynamics turns out to be periodic. As an 
example, Figs. 3(a) and (b) show a typical behavior of 
k„ and V n in an Hermitian lattice for parameter values 
N = 10, a = -N and 7 = 0, which yield E = M 2 - 21M 
according to Eq.(ll). Periodic quantum evolution is 
thus expected at the period T = ir. Self-imaging is 
demonstrated, as an example, in Fig. 3(c), where the 
behavior of the site occupation probability \ip(n = 0,t)| 2 
at the lattice edge n = versus time t for initial 
single-site excitation i/j(n, 0) = ^n,o is depicted. 



Conclusions. A novel class of truncated tight-binding 
lattices that sustain a periodic quantum evolution of the 
wave packet has been presented. In particular, we have 
shown the existence of Bloch oscillations in Hermitian 
truncated lattices, Bloch oscillations in non-Hermitian 
lattices with an imaginary site energy gradient, and 
discussed self-imaging phenomena in finite lattices with 
commensurate energy levels that depend quadratically 
on the quantum number. Photonic waveguide arrays, 
realized for instance by femtosecond laser writing in 
fused silica 0, could provide a possible experimental 
set-up to realize the Hermitian lattice models discussed 
in this paper (Figs.l and 3). In such arrayed structures, 
simultaneous engineering of coupling rates K n and site 
energies V n can be achieved by appropriate design of 
distances between adjacent waveguides (which determine 
the values of n n ) and choice of waveguide writing speed 
(which determines the values of V n ). 
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